3.5.55 \(\int \frac {\sqrt {a+b x^3} (c+d x+e x^2+f x^3+g x^4)}{x^7} \, dx\) [455]

3.5.55.1 Optimal result
3.5.55.2 Mathematica [C] (verified)
3.5.55.3 Rubi [A] (verified)
3.5.55.4 Maple [A] (verified)
3.5.55.5 Fricas [C] (verification not implemented)
3.5.55.6 Sympy [A] (verification not implemented)
3.5.55.7 Maxima [F]
3.5.55.8 Giac [F]
3.5.55.9 Mupad [F(-1)]

3.5.55.1 Optimal result

Integrand size = 35, antiderivative size = 659 \[ \int \frac {\sqrt {a+b x^3} \left (c+d x+e x^2+f x^3+g x^4\right )}{x^7} \, dx=-\frac {1}{60} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right ) \sqrt {a+b x^3}-\frac {b c \sqrt {a+b x^3}}{12 a x^3}-\frac {3 b d \sqrt {a+b x^3}}{20 a x^2}-\frac {3 b e \sqrt {a+b x^3}}{8 a x}+\frac {3 b^{4/3} e \sqrt {a+b x^3}}{8 a \left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )}+\frac {b (b c-4 a f) \text {arctanh}\left (\frac {\sqrt {a+b x^3}}{\sqrt {a}}\right )}{12 a^{3/2}}-\frac {3 \sqrt [4]{3} \sqrt {2-\sqrt {3}} b^{4/3} e \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} E\left (\arcsin \left (\frac {\left (1-\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}{\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}\right )|-7-4 \sqrt {3}\right )}{16 a^{2/3} \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}-\frac {3^{3/4} \sqrt {2+\sqrt {3}} b^{2/3} \left (2 b d+5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} b^{2/3} e-20 a g\right ) \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\left (1-\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}{\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x}\right ),-7-4 \sqrt {3}\right )}{40 a \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}} \]

output
1/12*b*(-4*a*f+b*c)*arctanh((b*x^3+a)^(1/2)/a^(1/2))/a^(3/2)-1/60*(10*c/x^ 
6+12*d/x^5+15*e/x^4+20*f/x^3+30*g/x^2)*(b*x^3+a)^(1/2)-1/12*b*c*(b*x^3+a)^ 
(1/2)/a/x^3-3/20*b*d*(b*x^3+a)^(1/2)/a/x^2-3/8*b*e*(b*x^3+a)^(1/2)/a/x+3/8 
*b^(4/3)*e*(b*x^3+a)^(1/2)/a/(b^(1/3)*x+a^(1/3)*(1+3^(1/2)))-3/16*3^(1/4)* 
b^(4/3)*e*(a^(1/3)+b^(1/3)*x)*EllipticE((b^(1/3)*x+a^(1/3)*(1-3^(1/2)))/(b 
^(1/3)*x+a^(1/3)*(1+3^(1/2))),I*3^(1/2)+2*I)*(1/2*6^(1/2)-1/2*2^(1/2))*((a 
^(2/3)-a^(1/3)*b^(1/3)*x+b^(2/3)*x^2)/(b^(1/3)*x+a^(1/3)*(1+3^(1/2)))^2)^( 
1/2)/a^(2/3)/(b*x^3+a)^(1/2)/(a^(1/3)*(a^(1/3)+b^(1/3)*x)/(b^(1/3)*x+a^(1/ 
3)*(1+3^(1/2)))^2)^(1/2)-1/40*3^(3/4)*b^(2/3)*(a^(1/3)+b^(1/3)*x)*Elliptic 
F((b^(1/3)*x+a^(1/3)*(1-3^(1/2)))/(b^(1/3)*x+a^(1/3)*(1+3^(1/2))),I*3^(1/2 
)+2*I)*(2*b*d-20*a*g+5*a^(1/3)*b^(2/3)*e*(1-3^(1/2)))*(1/2*6^(1/2)+1/2*2^( 
1/2))*((a^(2/3)-a^(1/3)*b^(1/3)*x+b^(2/3)*x^2)/(b^(1/3)*x+a^(1/3)*(1+3^(1/ 
2)))^2)^(1/2)/a/(b*x^3+a)^(1/2)/(a^(1/3)*(a^(1/3)+b^(1/3)*x)/(b^(1/3)*x+a^ 
(1/3)*(1+3^(1/2)))^2)^(1/2)
 
3.5.55.2 Mathematica [C] (verified)

Result contains higher order function than in optimal. Order 5 vs. order 4 in optimal.

Time = 10.43 (sec) , antiderivative size = 211, normalized size of antiderivative = 0.32 \[ \int \frac {\sqrt {a+b x^3} \left (c+d x+e x^2+f x^3+g x^4\right )}{x^7} \, dx=-\frac {\sqrt {a+b x^3} \left (36 a^3 d \operatorname {Hypergeometric2F1}\left (-\frac {5}{3},-\frac {1}{2},-\frac {2}{3},-\frac {b x^3}{a}\right )+5 x \left (9 a^3 e \operatorname {Hypergeometric2F1}\left (-\frac {4}{3},-\frac {1}{2},-\frac {1}{3},-\frac {b x^3}{a}\right )+2 x \left (6 a^2 f \left (a \sqrt {1+\frac {b x^3}{a}}+b x^3 \text {arctanh}\left (\sqrt {1+\frac {b x^3}{a}}\right )\right )+9 a^3 g x \operatorname {Hypergeometric2F1}\left (-\frac {2}{3},-\frac {1}{2},\frac {1}{3},-\frac {b x^3}{a}\right )+4 b^2 c x^3 \left (a+b x^3\right ) \sqrt {1+\frac {b x^3}{a}} \operatorname {Hypergeometric2F1}\left (\frac {3}{2},3,\frac {5}{2},1+\frac {b x^3}{a}\right )\right )\right )\right )}{180 a^3 x^5 \sqrt {1+\frac {b x^3}{a}}} \]

input
Integrate[(Sqrt[a + b*x^3]*(c + d*x + e*x^2 + f*x^3 + g*x^4))/x^7,x]
 
output
-1/180*(Sqrt[a + b*x^3]*(36*a^3*d*Hypergeometric2F1[-5/3, -1/2, -2/3, -((b 
*x^3)/a)] + 5*x*(9*a^3*e*Hypergeometric2F1[-4/3, -1/2, -1/3, -((b*x^3)/a)] 
 + 2*x*(6*a^2*f*(a*Sqrt[1 + (b*x^3)/a] + b*x^3*ArcTanh[Sqrt[1 + (b*x^3)/a] 
]) + 9*a^3*g*x*Hypergeometric2F1[-2/3, -1/2, 1/3, -((b*x^3)/a)] + 4*b^2*c* 
x^3*(a + b*x^3)*Sqrt[1 + (b*x^3)/a]*Hypergeometric2F1[3/2, 3, 5/2, 1 + (b* 
x^3)/a]))))/(a^3*x^5*Sqrt[1 + (b*x^3)/a])
 
3.5.55.3 Rubi [A] (verified)

Time = 1.29 (sec) , antiderivative size = 671, normalized size of antiderivative = 1.02, number of steps used = 15, number of rules used = 14, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.400, Rules used = {2364, 27, 2374, 27, 2374, 27, 2374, 2371, 798, 73, 221, 2417, 759, 2416}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\sqrt {a+b x^3} \left (c+d x+e x^2+f x^3+g x^4\right )}{x^7} \, dx\)

\(\Big \downarrow \) 2364

\(\displaystyle -\frac {3}{2} b \int -\frac {30 g x^4+20 f x^3+15 e x^2+12 d x+10 c}{60 x^4 \sqrt {b x^3+a}}dx-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {1}{40} b \int \frac {30 g x^4+20 f x^3+15 e x^2+12 d x+10 c}{x^4 \sqrt {b x^3+a}}dx-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 2374

\(\displaystyle \frac {1}{40} b \left (-\frac {\int -\frac {6 \left (30 a g x^3-5 (b c-4 a f) x^2+15 a e x+12 a d\right )}{x^3 \sqrt {b x^3+a}}dx}{6 a}-\frac {10 c \sqrt {a+b x^3}}{3 a x^3}\right )-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {1}{40} b \left (\frac {\int \frac {30 a g x^3-5 (b c-4 a f) x^2+15 a e x+12 a d}{x^3 \sqrt {b x^3+a}}dx}{a}-\frac {10 c \sqrt {a+b x^3}}{3 a x^3}\right )-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 2374

\(\displaystyle \frac {1}{40} b \left (\frac {-\frac {\int -\frac {4 \left (15 e a^2-3 (b d-10 a g) x^2 a-5 (b c-4 a f) x a\right )}{x^2 \sqrt {b x^3+a}}dx}{4 a}-\frac {6 d \sqrt {a+b x^3}}{x^2}}{a}-\frac {10 c \sqrt {a+b x^3}}{3 a x^3}\right )-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {1}{40} b \left (\frac {\frac {\int \frac {15 e a^2-3 (b d-10 a g) x^2 a-5 (b c-4 a f) x a}{x^2 \sqrt {b x^3+a}}dx}{a}-\frac {6 d \sqrt {a+b x^3}}{x^2}}{a}-\frac {10 c \sqrt {a+b x^3}}{3 a x^3}\right )-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 2374

\(\displaystyle \frac {1}{40} b \left (\frac {\frac {-\frac {\int \frac {-15 b e x^2 a^2+10 (b c-4 a f) a^2+6 (b d-10 a g) x a^2}{x \sqrt {b x^3+a}}dx}{2 a}-\frac {15 a e \sqrt {a+b x^3}}{x}}{a}-\frac {6 d \sqrt {a+b x^3}}{x^2}}{a}-\frac {10 c \sqrt {a+b x^3}}{3 a x^3}\right )-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 2371

\(\displaystyle \frac {1}{40} b \left (\frac {\frac {-\frac {10 a^2 (b c-4 a f) \int \frac {1}{x \sqrt {b x^3+a}}dx+\int \frac {6 a^2 (b d-10 a g)-15 a^2 b e x}{\sqrt {b x^3+a}}dx}{2 a}-\frac {15 a e \sqrt {a+b x^3}}{x}}{a}-\frac {6 d \sqrt {a+b x^3}}{x^2}}{a}-\frac {10 c \sqrt {a+b x^3}}{3 a x^3}\right )-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 798

\(\displaystyle \frac {1}{40} b \left (\frac {\frac {-\frac {\frac {10}{3} a^2 (b c-4 a f) \int \frac {1}{x^3 \sqrt {b x^3+a}}dx^3+\int \frac {6 a^2 (b d-10 a g)-15 a^2 b e x}{\sqrt {b x^3+a}}dx}{2 a}-\frac {15 a e \sqrt {a+b x^3}}{x}}{a}-\frac {6 d \sqrt {a+b x^3}}{x^2}}{a}-\frac {10 c \sqrt {a+b x^3}}{3 a x^3}\right )-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 73

\(\displaystyle \frac {1}{40} b \left (\frac {\frac {-\frac {\frac {20 a^2 (b c-4 a f) \int \frac {1}{\frac {x^6}{b}-\frac {a}{b}}d\sqrt {b x^3+a}}{3 b}+\int \frac {6 a^2 (b d-10 a g)-15 a^2 b e x}{\sqrt {b x^3+a}}dx}{2 a}-\frac {15 a e \sqrt {a+b x^3}}{x}}{a}-\frac {6 d \sqrt {a+b x^3}}{x^2}}{a}-\frac {10 c \sqrt {a+b x^3}}{3 a x^3}\right )-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 221

\(\displaystyle \frac {1}{40} b \left (\frac {\frac {-\frac {\int \frac {6 a^2 (b d-10 a g)-15 a^2 b e x}{\sqrt {b x^3+a}}dx-\frac {20}{3} a^{3/2} \text {arctanh}\left (\frac {\sqrt {a+b x^3}}{\sqrt {a}}\right ) (b c-4 a f)}{2 a}-\frac {15 a e \sqrt {a+b x^3}}{x}}{a}-\frac {6 d \sqrt {a+b x^3}}{x^2}}{a}-\frac {10 c \sqrt {a+b x^3}}{3 a x^3}\right )-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 2417

\(\displaystyle \frac {1}{40} b \left (\frac {\frac {-\frac {3 a^2 \left (5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} b^{2/3} e-20 a g+2 b d\right ) \int \frac {1}{\sqrt {b x^3+a}}dx-15 a^2 b^{2/3} e \int \frac {\sqrt [3]{b} x+\left (1-\sqrt {3}\right ) \sqrt [3]{a}}{\sqrt {b x^3+a}}dx-\frac {20}{3} a^{3/2} \text {arctanh}\left (\frac {\sqrt {a+b x^3}}{\sqrt {a}}\right ) (b c-4 a f)}{2 a}-\frac {15 a e \sqrt {a+b x^3}}{x}}{a}-\frac {6 d \sqrt {a+b x^3}}{x^2}}{a}-\frac {10 c \sqrt {a+b x^3}}{3 a x^3}\right )-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 759

\(\displaystyle \frac {1}{40} b \left (\frac {\frac {-\frac {-15 a^2 b^{2/3} e \int \frac {\sqrt [3]{b} x+\left (1-\sqrt {3}\right ) \sqrt [3]{a}}{\sqrt {b x^3+a}}dx-\frac {20}{3} a^{3/2} \text {arctanh}\left (\frac {\sqrt {a+b x^3}}{\sqrt {a}}\right ) (b c-4 a f)+\frac {2\ 3^{3/4} \sqrt {2+\sqrt {3}} a^2 \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt [3]{b} x+\left (1-\sqrt {3}\right ) \sqrt [3]{a}}{\sqrt [3]{b} x+\left (1+\sqrt {3}\right ) \sqrt [3]{a}}\right ),-7-4 \sqrt {3}\right ) \left (5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} b^{2/3} e-20 a g+2 b d\right )}{\sqrt [3]{b} \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}}{2 a}-\frac {15 a e \sqrt {a+b x^3}}{x}}{a}-\frac {6 d \sqrt {a+b x^3}}{x^2}}{a}-\frac {10 c \sqrt {a+b x^3}}{3 a x^3}\right )-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

\(\Big \downarrow \) 2416

\(\displaystyle \frac {1}{40} b \left (\frac {\frac {-\frac {-\frac {20}{3} a^{3/2} \text {arctanh}\left (\frac {\sqrt {a+b x^3}}{\sqrt {a}}\right ) (b c-4 a f)+\frac {2\ 3^{3/4} \sqrt {2+\sqrt {3}} a^2 \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt [3]{b} x+\left (1-\sqrt {3}\right ) \sqrt [3]{a}}{\sqrt [3]{b} x+\left (1+\sqrt {3}\right ) \sqrt [3]{a}}\right ),-7-4 \sqrt {3}\right ) \left (5 \left (1-\sqrt {3}\right ) \sqrt [3]{a} b^{2/3} e-20 a g+2 b d\right )}{\sqrt [3]{b} \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}-15 a^2 b^{2/3} e \left (\frac {2 \sqrt {a+b x^3}}{\sqrt [3]{b} \left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )}-\frac {\sqrt [4]{3} \sqrt {2-\sqrt {3}} \sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right ) \sqrt {\frac {a^{2/3}-\sqrt [3]{a} \sqrt [3]{b} x+b^{2/3} x^2}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} E\left (\arcsin \left (\frac {\sqrt [3]{b} x+\left (1-\sqrt {3}\right ) \sqrt [3]{a}}{\sqrt [3]{b} x+\left (1+\sqrt {3}\right ) \sqrt [3]{a}}\right )|-7-4 \sqrt {3}\right )}{\sqrt [3]{b} \sqrt {\frac {\sqrt [3]{a} \left (\sqrt [3]{a}+\sqrt [3]{b} x\right )}{\left (\left (1+\sqrt {3}\right ) \sqrt [3]{a}+\sqrt [3]{b} x\right )^2}} \sqrt {a+b x^3}}\right )}{2 a}-\frac {15 a e \sqrt {a+b x^3}}{x}}{a}-\frac {6 d \sqrt {a+b x^3}}{x^2}}{a}-\frac {10 c \sqrt {a+b x^3}}{3 a x^3}\right )-\frac {1}{60} \sqrt {a+b x^3} \left (\frac {10 c}{x^6}+\frac {12 d}{x^5}+\frac {15 e}{x^4}+\frac {20 f}{x^3}+\frac {30 g}{x^2}\right )\)

input
Int[(Sqrt[a + b*x^3]*(c + d*x + e*x^2 + f*x^3 + g*x^4))/x^7,x]
 
output
-1/60*(((10*c)/x^6 + (12*d)/x^5 + (15*e)/x^4 + (20*f)/x^3 + (30*g)/x^2)*Sq 
rt[a + b*x^3]) + (b*((-10*c*Sqrt[a + b*x^3])/(3*a*x^3) + ((-6*d*Sqrt[a + b 
*x^3])/x^2 + ((-15*a*e*Sqrt[a + b*x^3])/x - ((-20*a^(3/2)*(b*c - 4*a*f)*Ar 
cTanh[Sqrt[a + b*x^3]/Sqrt[a]])/3 - 15*a^2*b^(2/3)*e*((2*Sqrt[a + b*x^3])/ 
(b^(1/3)*((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)) - (3^(1/4)*Sqrt[2 - Sqrt[3]] 
*a^(1/3)*(a^(1/3) + b^(1/3)*x)*Sqrt[(a^(2/3) - a^(1/3)*b^(1/3)*x + b^(2/3) 
*x^2)/((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)^2]*EllipticE[ArcSin[((1 - Sqrt[3 
])*a^(1/3) + b^(1/3)*x)/((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)], -7 - 4*Sqrt[ 
3]])/(b^(1/3)*Sqrt[(a^(1/3)*(a^(1/3) + b^(1/3)*x))/((1 + Sqrt[3])*a^(1/3) 
+ b^(1/3)*x)^2]*Sqrt[a + b*x^3])) + (2*3^(3/4)*Sqrt[2 + Sqrt[3]]*a^2*(2*b* 
d + 5*(1 - Sqrt[3])*a^(1/3)*b^(2/3)*e - 20*a*g)*(a^(1/3) + b^(1/3)*x)*Sqrt 
[(a^(2/3) - a^(1/3)*b^(1/3)*x + b^(2/3)*x^2)/((1 + Sqrt[3])*a^(1/3) + b^(1 
/3)*x)^2]*EllipticF[ArcSin[((1 - Sqrt[3])*a^(1/3) + b^(1/3)*x)/((1 + Sqrt[ 
3])*a^(1/3) + b^(1/3)*x)], -7 - 4*Sqrt[3]])/(b^(1/3)*Sqrt[(a^(1/3)*(a^(1/3 
) + b^(1/3)*x))/((1 + Sqrt[3])*a^(1/3) + b^(1/3)*x)^2]*Sqrt[a + b*x^3]))/( 
2*a))/a)/a))/40
 

3.5.55.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 73
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> With[ 
{p = Denominator[m]}, Simp[p/b   Subst[Int[x^(p*(m + 1) - 1)*(c - a*(d/b) + 
 d*(x^p/b))^n, x], x, (a + b*x)^(1/p)], x]] /; FreeQ[{a, b, c, d}, x] && Lt 
Q[-1, m, 0] && LeQ[-1, n, 0] && LeQ[Denominator[n], Denominator[m]] && IntL 
inearQ[a, b, c, d, m, n, x]
 

rule 221
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x 
/Rt[-a/b, 2]], x] /; FreeQ[{a, b}, x] && NegQ[a/b]
 

rule 759
Int[1/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = Numer[Rt[b/a, 3]], 
s = Denom[Rt[b/a, 3]]}, Simp[2*Sqrt[2 + Sqrt[3]]*(s + r*x)*(Sqrt[(s^2 - r*s 
*x + r^2*x^2)/((1 + Sqrt[3])*s + r*x)^2]/(3^(1/4)*r*Sqrt[a + b*x^3]*Sqrt[s* 
((s + r*x)/((1 + Sqrt[3])*s + r*x)^2)]))*EllipticF[ArcSin[((1 - Sqrt[3])*s 
+ r*x)/((1 + Sqrt[3])*s + r*x)], -7 - 4*Sqrt[3]], x]] /; FreeQ[{a, b}, x] & 
& PosQ[a]
 

rule 798
Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[1/n   Subst 
[Int[x^(Simplify[(m + 1)/n] - 1)*(a + b*x)^p, x], x, x^n], x] /; FreeQ[{a, 
b, m, n, p}, x] && IntegerQ[Simplify[(m + 1)/n]]
 

rule 2364
Int[(Pq_)*(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_.))^(p_), x_Symbol] :> Module[{u 
 = IntHide[x^m*Pq, x]}, Simp[u*(a + b*x^n)^p, x] - Simp[b*n*p   Int[x^(m + 
n)*(a + b*x^n)^(p - 1)*ExpandToSum[u/x^(m + 1), x], x], x]] /; FreeQ[{a, b} 
, x] && PolyQ[Pq, x] && IGtQ[n, 0] && GtQ[p, 0] && LtQ[m + Expon[Pq, x] + 1 
, 0]
 

rule 2371
Int[(Pq_)/((x_)*Sqrt[(a_) + (b_.)*(x_)^(n_)]), x_Symbol] :> Simp[Coeff[Pq, 
x, 0]   Int[1/(x*Sqrt[a + b*x^n]), x], x] + Int[ExpandToSum[(Pq - Coeff[Pq, 
 x, 0])/x, x]/Sqrt[a + b*x^n], x] /; FreeQ[{a, b}, x] && PolyQ[Pq, x] && IG 
tQ[n, 0] && NeQ[Coeff[Pq, x, 0], 0]
 

rule 2374
Int[(Pq_)*((c_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Wit 
h[{Pq0 = Coeff[Pq, x, 0]}, Simp[Pq0*(c*x)^(m + 1)*((a + b*x^n)^(p + 1)/(a*c 
*(m + 1))), x] + Simp[1/(2*a*c*(m + 1))   Int[(c*x)^(m + 1)*ExpandToSum[2*a 
*(m + 1)*((Pq - Pq0)/x) - 2*b*Pq0*(m + n*(p + 1) + 1)*x^(n - 1), x]*(a + b* 
x^n)^p, x], x] /; NeQ[Pq0, 0]] /; FreeQ[{a, b, c, p}, x] && PolyQ[Pq, x] && 
 IGtQ[n, 0] && LtQ[m, -1] && LeQ[n - 1, Expon[Pq, x]]
 

rule 2416
Int[((c_) + (d_.)*(x_))/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = N 
umer[Simplify[(1 - Sqrt[3])*(d/c)]], s = Denom[Simplify[(1 - Sqrt[3])*(d/c) 
]]}, Simp[2*d*s^3*(Sqrt[a + b*x^3]/(a*r^2*((1 + Sqrt[3])*s + r*x))), x] - S 
imp[3^(1/4)*Sqrt[2 - Sqrt[3]]*d*s*(s + r*x)*(Sqrt[(s^2 - r*s*x + r^2*x^2)/( 
(1 + Sqrt[3])*s + r*x)^2]/(r^2*Sqrt[a + b*x^3]*Sqrt[s*((s + r*x)/((1 + Sqrt 
[3])*s + r*x)^2)]))*EllipticE[ArcSin[((1 - Sqrt[3])*s + r*x)/((1 + Sqrt[3]) 
*s + r*x)], -7 - 4*Sqrt[3]], x]] /; FreeQ[{a, b, c, d}, x] && PosQ[a] && Eq 
Q[b*c^3 - 2*(5 - 3*Sqrt[3])*a*d^3, 0]
 

rule 2417
Int[((c_) + (d_.)*(x_))/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = N 
umer[Rt[b/a, 3]], s = Denom[Rt[b/a, 3]]}, Simp[(c*r - (1 - Sqrt[3])*d*s)/r 
  Int[1/Sqrt[a + b*x^3], x], x] + Simp[d/r   Int[((1 - Sqrt[3])*s + r*x)/Sq 
rt[a + b*x^3], x], x]] /; FreeQ[{a, b, c, d}, x] && PosQ[a] && NeQ[b*c^3 - 
2*(5 - 3*Sqrt[3])*a*d^3, 0]
 
3.5.55.4 Maple [A] (verified)

Time = 1.82 (sec) , antiderivative size = 883, normalized size of antiderivative = 1.34

method result size
elliptic \(\text {Expression too large to display}\) \(883\)
risch \(\text {Expression too large to display}\) \(1102\)
default \(\text {Expression too large to display}\) \(1180\)

input
int((g*x^4+f*x^3+e*x^2+d*x+c)*(b*x^3+a)^(1/2)/x^7,x,method=_RETURNVERBOSE)
 
output
-1/6*c*(b*x^3+a)^(1/2)/x^6-1/5*d*(b*x^3+a)^(1/2)/x^5-1/4*e*(b*x^3+a)^(1/2) 
/x^4-1/12*(4*a*f+b*c)/a*(b*x^3+a)^(1/2)/x^3-1/20/a*(10*a*g+3*b*d)*(b*x^3+a 
)^(1/2)/x^2-3/8*b*e*(b*x^3+a)^(1/2)/a/x-2/3*I*(g*b-1/40*b/a*(10*a*g+3*b*d) 
)*3^(1/2)/b*(-a*b^2)^(1/3)*(I*(x+1/2/b*(-a*b^2)^(1/3)-1/2*I*3^(1/2)/b*(-a* 
b^2)^(1/3))*3^(1/2)*b/(-a*b^2)^(1/3))^(1/2)*((x-1/b*(-a*b^2)^(1/3))/(-3/2/ 
b*(-a*b^2)^(1/3)+1/2*I*3^(1/2)/b*(-a*b^2)^(1/3)))^(1/2)*(-I*(x+1/2/b*(-a*b 
^2)^(1/3)+1/2*I*3^(1/2)/b*(-a*b^2)^(1/3))*3^(1/2)*b/(-a*b^2)^(1/3))^(1/2)/ 
(b*x^3+a)^(1/2)*EllipticF(1/3*3^(1/2)*(I*(x+1/2/b*(-a*b^2)^(1/3)-1/2*I*3^( 
1/2)/b*(-a*b^2)^(1/3))*3^(1/2)*b/(-a*b^2)^(1/3))^(1/2),(I*3^(1/2)/b*(-a*b^ 
2)^(1/3)/(-3/2/b*(-a*b^2)^(1/3)+1/2*I*3^(1/2)/b*(-a*b^2)^(1/3)))^(1/2))-1/ 
8*I/a*b*e*3^(1/2)*(-a*b^2)^(1/3)*(I*(x+1/2/b*(-a*b^2)^(1/3)-1/2*I*3^(1/2)/ 
b*(-a*b^2)^(1/3))*3^(1/2)*b/(-a*b^2)^(1/3))^(1/2)*((x-1/b*(-a*b^2)^(1/3))/ 
(-3/2/b*(-a*b^2)^(1/3)+1/2*I*3^(1/2)/b*(-a*b^2)^(1/3)))^(1/2)*(-I*(x+1/2/b 
*(-a*b^2)^(1/3)+1/2*I*3^(1/2)/b*(-a*b^2)^(1/3))*3^(1/2)*b/(-a*b^2)^(1/3))^ 
(1/2)/(b*x^3+a)^(1/2)*((-3/2/b*(-a*b^2)^(1/3)+1/2*I*3^(1/2)/b*(-a*b^2)^(1/ 
3))*EllipticE(1/3*3^(1/2)*(I*(x+1/2/b*(-a*b^2)^(1/3)-1/2*I*3^(1/2)/b*(-a*b 
^2)^(1/3))*3^(1/2)*b/(-a*b^2)^(1/3))^(1/2),(I*3^(1/2)/b*(-a*b^2)^(1/3)/(-3 
/2/b*(-a*b^2)^(1/3)+1/2*I*3^(1/2)/b*(-a*b^2)^(1/3)))^(1/2))+1/b*(-a*b^2)^( 
1/3)*EllipticF(1/3*3^(1/2)*(I*(x+1/2/b*(-a*b^2)^(1/3)-1/2*I*3^(1/2)/b*(-a* 
b^2)^(1/3))*3^(1/2)*b/(-a*b^2)^(1/3))^(1/2),(I*3^(1/2)/b*(-a*b^2)^(1/3)...
 
3.5.55.5 Fricas [C] (verification not implemented)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 0.17 (sec) , antiderivative size = 404, normalized size of antiderivative = 0.61 \[ \int \frac {\sqrt {a+b x^3} \left (c+d x+e x^2+f x^3+g x^4\right )}{x^7} \, dx=\left [-\frac {90 \, a b^{\frac {3}{2}} e x^{6} {\rm weierstrassZeta}\left (0, -\frac {4 \, a}{b}, {\rm weierstrassPInverse}\left (0, -\frac {4 \, a}{b}, x\right )\right ) + 5 \, {\left (b^{2} c - 4 \, a b f\right )} \sqrt {a} x^{6} \log \left (\frac {b^{2} x^{6} + 8 \, a b x^{3} - 4 \, {\left (b x^{3} + 2 \, a\right )} \sqrt {b x^{3} + a} \sqrt {a} + 8 \, a^{2}}{x^{6}}\right ) + 36 \, {\left (a b d - 10 \, a^{2} g\right )} \sqrt {b} x^{6} {\rm weierstrassPInverse}\left (0, -\frac {4 \, a}{b}, x\right ) + 2 \, {\left (45 \, a b e x^{5} + 30 \, a^{2} e x^{2} + 6 \, {\left (3 \, a b d + 10 \, a^{2} g\right )} x^{4} + 24 \, a^{2} d x + 10 \, {\left (a b c + 4 \, a^{2} f\right )} x^{3} + 20 \, a^{2} c\right )} \sqrt {b x^{3} + a}}{240 \, a^{2} x^{6}}, -\frac {45 \, a b^{\frac {3}{2}} e x^{6} {\rm weierstrassZeta}\left (0, -\frac {4 \, a}{b}, {\rm weierstrassPInverse}\left (0, -\frac {4 \, a}{b}, x\right )\right ) + 5 \, {\left (b^{2} c - 4 \, a b f\right )} \sqrt {-a} x^{6} \arctan \left (\frac {{\left (b x^{3} + 2 \, a\right )} \sqrt {b x^{3} + a} \sqrt {-a}}{2 \, {\left (a b x^{3} + a^{2}\right )}}\right ) + 18 \, {\left (a b d - 10 \, a^{2} g\right )} \sqrt {b} x^{6} {\rm weierstrassPInverse}\left (0, -\frac {4 \, a}{b}, x\right ) + {\left (45 \, a b e x^{5} + 30 \, a^{2} e x^{2} + 6 \, {\left (3 \, a b d + 10 \, a^{2} g\right )} x^{4} + 24 \, a^{2} d x + 10 \, {\left (a b c + 4 \, a^{2} f\right )} x^{3} + 20 \, a^{2} c\right )} \sqrt {b x^{3} + a}}{120 \, a^{2} x^{6}}\right ] \]

input
integrate((g*x^4+f*x^3+e*x^2+d*x+c)*(b*x^3+a)^(1/2)/x^7,x, algorithm="fric 
as")
 
output
[-1/240*(90*a*b^(3/2)*e*x^6*weierstrassZeta(0, -4*a/b, weierstrassPInverse 
(0, -4*a/b, x)) + 5*(b^2*c - 4*a*b*f)*sqrt(a)*x^6*log((b^2*x^6 + 8*a*b*x^3 
 - 4*(b*x^3 + 2*a)*sqrt(b*x^3 + a)*sqrt(a) + 8*a^2)/x^6) + 36*(a*b*d - 10* 
a^2*g)*sqrt(b)*x^6*weierstrassPInverse(0, -4*a/b, x) + 2*(45*a*b*e*x^5 + 3 
0*a^2*e*x^2 + 6*(3*a*b*d + 10*a^2*g)*x^4 + 24*a^2*d*x + 10*(a*b*c + 4*a^2* 
f)*x^3 + 20*a^2*c)*sqrt(b*x^3 + a))/(a^2*x^6), -1/120*(45*a*b^(3/2)*e*x^6* 
weierstrassZeta(0, -4*a/b, weierstrassPInverse(0, -4*a/b, x)) + 5*(b^2*c - 
 4*a*b*f)*sqrt(-a)*x^6*arctan(1/2*(b*x^3 + 2*a)*sqrt(b*x^3 + a)*sqrt(-a)/( 
a*b*x^3 + a^2)) + 18*(a*b*d - 10*a^2*g)*sqrt(b)*x^6*weierstrassPInverse(0, 
 -4*a/b, x) + (45*a*b*e*x^5 + 30*a^2*e*x^2 + 6*(3*a*b*d + 10*a^2*g)*x^4 + 
24*a^2*d*x + 10*(a*b*c + 4*a^2*f)*x^3 + 20*a^2*c)*sqrt(b*x^3 + a))/(a^2*x^ 
6)]
 
3.5.55.6 Sympy [A] (verification not implemented)

Time = 4.76 (sec) , antiderivative size = 304, normalized size of antiderivative = 0.46 \[ \int \frac {\sqrt {a+b x^3} \left (c+d x+e x^2+f x^3+g x^4\right )}{x^7} \, dx=\frac {\sqrt {a} d \Gamma \left (- \frac {5}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {5}{3}, - \frac {1}{2} \\ - \frac {2}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 x^{5} \Gamma \left (- \frac {2}{3}\right )} + \frac {\sqrt {a} e \Gamma \left (- \frac {4}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {4}{3}, - \frac {1}{2} \\ - \frac {1}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 x^{4} \Gamma \left (- \frac {1}{3}\right )} + \frac {\sqrt {a} g \Gamma \left (- \frac {2}{3}\right ) {{}_{2}F_{1}\left (\begin {matrix} - \frac {2}{3}, - \frac {1}{2} \\ \frac {1}{3} \end {matrix}\middle | {\frac {b x^{3} e^{i \pi }}{a}} \right )}}{3 x^{2} \Gamma \left (\frac {1}{3}\right )} - \frac {a c}{6 \sqrt {b} x^{\frac {15}{2}} \sqrt {\frac {a}{b x^{3}} + 1}} - \frac {\sqrt {b} c}{4 x^{\frac {9}{2}} \sqrt {\frac {a}{b x^{3}} + 1}} - \frac {\sqrt {b} f \sqrt {\frac {a}{b x^{3}} + 1}}{3 x^{\frac {3}{2}}} - \frac {b^{\frac {3}{2}} c}{12 a x^{\frac {3}{2}} \sqrt {\frac {a}{b x^{3}} + 1}} - \frac {b f \operatorname {asinh}{\left (\frac {\sqrt {a}}{\sqrt {b} x^{\frac {3}{2}}} \right )}}{3 \sqrt {a}} + \frac {b^{2} c \operatorname {asinh}{\left (\frac {\sqrt {a}}{\sqrt {b} x^{\frac {3}{2}}} \right )}}{12 a^{\frac {3}{2}}} \]

input
integrate((g*x**4+f*x**3+e*x**2+d*x+c)*(b*x**3+a)**(1/2)/x**7,x)
 
output
sqrt(a)*d*gamma(-5/3)*hyper((-5/3, -1/2), (-2/3,), b*x**3*exp_polar(I*pi)/ 
a)/(3*x**5*gamma(-2/3)) + sqrt(a)*e*gamma(-4/3)*hyper((-4/3, -1/2), (-1/3, 
), b*x**3*exp_polar(I*pi)/a)/(3*x**4*gamma(-1/3)) + sqrt(a)*g*gamma(-2/3)* 
hyper((-2/3, -1/2), (1/3,), b*x**3*exp_polar(I*pi)/a)/(3*x**2*gamma(1/3)) 
- a*c/(6*sqrt(b)*x**(15/2)*sqrt(a/(b*x**3) + 1)) - sqrt(b)*c/(4*x**(9/2)*s 
qrt(a/(b*x**3) + 1)) - sqrt(b)*f*sqrt(a/(b*x**3) + 1)/(3*x**(3/2)) - b**(3 
/2)*c/(12*a*x**(3/2)*sqrt(a/(b*x**3) + 1)) - b*f*asinh(sqrt(a)/(sqrt(b)*x* 
*(3/2)))/(3*sqrt(a)) + b**2*c*asinh(sqrt(a)/(sqrt(b)*x**(3/2)))/(12*a**(3/ 
2))
 
3.5.55.7 Maxima [F]

\[ \int \frac {\sqrt {a+b x^3} \left (c+d x+e x^2+f x^3+g x^4\right )}{x^7} \, dx=\int { \frac {{\left (g x^{4} + f x^{3} + e x^{2} + d x + c\right )} \sqrt {b x^{3} + a}}{x^{7}} \,d x } \]

input
integrate((g*x^4+f*x^3+e*x^2+d*x+c)*(b*x^3+a)^(1/2)/x^7,x, algorithm="maxi 
ma")
 
output
-1/24*(b^2*log((sqrt(b*x^3 + a) - sqrt(a))/(sqrt(b*x^3 + a) + sqrt(a)))/a^ 
(3/2) + 2*((b*x^3 + a)^(3/2)*b^2 + sqrt(b*x^3 + a)*a*b^2)/((b*x^3 + a)^2*a 
 - 2*(b*x^3 + a)*a^2 + a^3))*c + integrate(sqrt(b*x^3 + a)*(g*x^3 + f*x^2 
+ e*x + d)/x^6, x)
 
3.5.55.8 Giac [F]

\[ \int \frac {\sqrt {a+b x^3} \left (c+d x+e x^2+f x^3+g x^4\right )}{x^7} \, dx=\int { \frac {{\left (g x^{4} + f x^{3} + e x^{2} + d x + c\right )} \sqrt {b x^{3} + a}}{x^{7}} \,d x } \]

input
integrate((g*x^4+f*x^3+e*x^2+d*x+c)*(b*x^3+a)^(1/2)/x^7,x, algorithm="giac 
")
 
output
integrate((g*x^4 + f*x^3 + e*x^2 + d*x + c)*sqrt(b*x^3 + a)/x^7, x)
 
3.5.55.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\sqrt {a+b x^3} \left (c+d x+e x^2+f x^3+g x^4\right )}{x^7} \, dx=\int \frac {\sqrt {b\,x^3+a}\,\left (g\,x^4+f\,x^3+e\,x^2+d\,x+c\right )}{x^7} \,d x \]

input
int(((a + b*x^3)^(1/2)*(c + d*x + e*x^2 + f*x^3 + g*x^4))/x^7,x)
 
output
int(((a + b*x^3)^(1/2)*(c + d*x + e*x^2 + f*x^3 + g*x^4))/x^7, x)